Exact solution of the open Heisenberg chain with two impurities 
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We propose an integrable model of the spin- 1/2 Heisenberg chain coupled to two impurity mo- 
. . . ments. With the open boundary conditions at the impurity sites, the model can be exactly solved 

00 ■ for arbitrary impurity spin and arbitrary exchange constants between the bulk and the impurities. 

The absence of redundant terms in the hamiltonian makes the model very reasonable. The hamil- 
tonian is diagonalized via algebraic Bethe ansatz. It is found that the impurity spins can only be 
screened (partially for S > 1/2) for antiferromagnetic coupling between the impurity and the bulk. 
Otherwise the impurity spins can not be screened. The residual entropy of the ground state and the 
Kondo temperature are also derived explicitly based on the thermodynamic Bethe ansatz and the 
local Fermi liquid theory. 

O ■ 75.30.Hx, 75.10.Jm, 05.50.+q 
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^ ■ I. INTRODUCTION 

, Recently, considerable attention has been focused on the problem of impurities embedded in a quantum chain. Using 
simpler bosonization and renormalization group techniques, Kane and Fisher have shown that a potential scattering 
center in a Luttinger liquid is driven to a strong-coupling fixed point by the repulsive interaction in the bulktl. This 

■ is the first time to show that a single impurity in a one-dimensional quantum system behaves rather different from 
' that in a Fermi liquid, and directly stimulates the study oitjihe problem of local perturbations to a Luttinger liquid 

(— I and especially on the Kondo problem in a Luttinger liquidHiS. It is jwell known now that the spin dynamics of the 
Q Kondo problem is equivalent to that of a spin chain with an impurityQ. 

O The integrable impuriigj- problem of the Heisenberg chain with periodic boundary condition was first considered by 

Andrei and JohannessonQ. They studied the integrable case of a spin S > 1/2 embedded in a spin-1/2 Heisenberg 
I chain. Subsequently, the problem was generalized to arbitrary spins by Lee and Schlottmann, and SchlottmaanEl. 
^ I Now the quantum inverse scattering method becomes a standard method to construcLintegrable impurity models. A 
recent example is the integrable impurity problem in the supersymmetric t — J mode]|13. In the QISM, the hamiltonian 
of the model is usually written as the logarithmic derivative of an homogeneous transfer matrix at a special value 
of the spectral parameter. Inputting some inhomogeneous vertex matrices (provided these matrices satisfy the same 
, Yang-Baxter relation of the homogeneous ones), then we obtain an inhomogeneous transfer matrix. Its logarithmic 
QQ derivative at some special value of the spectral parameter thus gives a hamiltonian with local interactions. However, 
to construct an integrable impurity model with periodic boundary condition, there is a prize to pay, namely some 
unphysical terms must present in the hamiltonian, though they may be irrelevant. Their presence is required by the 
intcgrablility. 

^ ,. In another hand, the open boundary problem for the quantum chain has been renewed due to Kane and Fisher's 
I ' observationEJ. It is found that the open boundary theory is vew |Useful to formulate both the thermodynamics and 
^ , the transport properties of the quantipa chains with impurityllIHl3. The open boundary problem of the integrable 
^ models was first considered by GaudinU, who studied the nonlinear Schrodinger model and the spin-1/2 Heiseaberg 
chain with simple open boundaries. Subsequently, hia_oiethod was generalized to the Hubbard model by Schulzt^ and 
the spin chain with boundary fields by Alcaraz et alJij. In addition, Sklyannin formulated the saHie result based on 
the QISM. In his theory, a new relation which now is called as the reflecting Yang-Baxter relationO was used. 

In this paper, we study the problem of an open spinl/2 Heisenberg chain coupled with two impurity spins sited at 

■ the ends of the system. The hamiltonian we shall consider reads 

H — - J ^ (Tj ■ (jj+i + Jr(jn ■ Sr + Jl<Ji ■ Sl, (1) 

where (Tj are the Pauli matrices; Sr^l are the impurity moments with an arbitrary spin S] N is the site number of 
the bulk; J is a positive constant and Jr^l are two arbitrary real constants which describe the coupling between the 
bulk and the impurities. We remark that the present model is very reasonable for the absence of redundant terms 
which can not be justified on physical grounds. 
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The structure of the present paper is the following: In the subsequent section, we construct the transfer matrix 
corresponding to the hamiltonian (1), thus the integrability of the present model can be directly justified. Based on 
the QISM, the Bethe ansatz equation and the eigenvalue of the hamiltonian will be derived. In sect. Ill, the ground 
state properties for different regions of parameter Jr,l will be discussed. Sect. IV is attributed to the residual entropy 
of the ground state and the low temperature specific heat. Concluding remarks will be given in sect.^. 



II. ALGEBRAIC BETHE ANSATZ 



In the frameworkjpf QISM, the integrable hamiltonian with open boundary condition is usually obtained from a 
monodromy matrixll3 

C/(A) - K^i\)T{\)K+i\)T-\-X), (2) 
where K±(X) are the reflecting matrices which satisfy the reflecting Yang-Baxter relation 

Si2{X - ^l)Ki{X)Sl2{X + /^)^±(Ai) = Kl{^i)Si2{X + ^l)Ki(X)S,2{X - m), (3) 

12 

with and acting on the spiace Vi,2 and Vi V2 respectively. 5i2 is the scattering matrix which satisfy 

the traditional Yang-Baxter relationll3. T(A) is the monodromy matrix for the periodic system which satisfy the 
Yang-Baxter relation 

Si2{X - m)Ti(A)T2(m) = ^2(^)7^1 (A) 5i2 (A - m). (4) 

As demonstrated by Sklyannin0, U{X) satisfies the same refiecting Yang-Baxter relation (3) as K± do. Thus the trace 
of U{X) gives an infinite number of conserved quantities. We remark that Sklyannin and many other authors used 
c- number K±{X) to construct their models, whjsre K±{X) only induce the boundary fields. The operator ones were 
first used by the present author and coworkersu to study the Kondo problem in one-dimensional strongly correlated 
electron systems. 

To construct the algebraic Bethe ansatz of the present model, we define 

TAX) = L+AX)LNriX)LN-lriX) ■ ■ ■ LiriX), 

fr{X) = LiriX)L2r{X) ■ ■ ■ Xat, ( A)i^, ( A) , (5) 

where Lji-(A) = iA + 1/2(1-1- ct-, • f), andi^^(A) = zA ± -|- | -|- r- S'i? with r an auxiliary Pauli matrix. Furthermore, 
we put 

i^-(A) = 1, 

A-+(A) - {tX + CL + ^+f-SLm-CL + l+r-SL), (6) 

where cr^l are two constants. Obviously, K±{X) defined in (6) satisfy the reflecting relation (3). The monodromy 
matrix U{X) is defined for our model as 

UriX) = A'_(A)r,(A)i^+(A)f,(A), (7) 

which satisfy the reflecting equation 

Sl2{X - ^)Ur,WSi2iX + M)C/r.(/i) = UMS^X + fl)Ur,WSi2{X ~ fl) , (8) 

with 5'i2(A ± fi) = Lr-^r2{X ± /i). The hamiltonian (1) is obtained by the following relation 

-iJ dX{X) TV 1 J 



H — 77^ — ^ 7^ Ia=0 —J — 



2 E^ is+hY'^cr 



4n.=K,L[(5 + i)2-c2] dX ^^-^ 2 2 ^^^(5^^, 

where X{X) ~ trrUr{X) and Jr^l are parametrized by cr^l as Jr^l — J/[[S + 1/2)^ ~ c^r l\- 

Although the model can be solved for arbitrary J^.l, we consider only the cr = cl ^ c (Jr = Jl = Ji) case in 
this paper. The general case can be formulated following the same procedure without any difficulty. We introduce 
the notation 
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U(X)-( \ 

. Some useful commutation relation can be formulated from (8) as 

[A{X),A{^)] = [A(A),Z)(m)] = [BiX),Bi^i)]=0, 



B{X)Ai^^), 



{X + n-i){2fi-i) 
2i{X- 
(A-/z)(2a 
4i(A - 



-i?(A)5(M), (10) 



{2ii-i){X + ii-i) 

where D{X) is defined as -D(A) = (2A — i)D{X) + iA{X). Therefore, the trace of Ur{X) can be expressed as 

X{X) = TrrUriX) = ^5(A) + ^^^(A)- (11) 
Define the pseudo vacuum state |0 > as 

cr+ 10 >= 5+ 10 >= 0. (12) 
The elements of Ur{X) acting on the pseudo vacuum state behave as 

C(A)|0 >=0, 

A{X)\0 >= a(A)|0 >= [{iX + S+^f- c']''{iX + lf^\0 >, (13) 

£'(A)|0 >= J(A)|0 >= {iXf^[{iX -S+^f- c2]2|0 > . 

Therefore, X{X) can be treated as the generating functional of an infinite number of conserved quantities (including 
the hamiltonian) and -B(A) is the creation operator of their eigenstates. An eigenstate of X{X) with M spins down 
can be constructed as 

M 

\Q>=l[B{Xj)\0>. (14) 

With the relations (10) and (13), we obtain the eigenvalue of X{X) acting on the state \Cl > as 

X(A)|0>=A(A;Ai,---,Am)|0>, 

^(x■x X ^- ^^^^ ^ (a-a, -z)(a + a,--2z) 

A(A, A„ • • • , Am) - ^ U (A-A,)(A + A,-z) ^''^ 

2(A - i)a(A) -f^ (A + Aj)(A - Aj + i) 
+ 2A-i U (A + A, -i)(A-A,)- 

However, the spectral parameters Xj are not independent each other but satisiy the following Bethe ansatz equation 

l^j + '^\'2N ^Xj — ic + iS .2 ^ Xj + ic + iS .2 _ -pr T-j-Aj— rAj+i 

^A~^^ ^Xj-ic-iS> ^Xj+ic-iS' - ^ll^llxj-rXt-i- ^ > 

The eigenvalue of the hamiltonian (1) acting on the state \fl > can be obtained from (9) and (15) as 

M 

^(Ai, ■ . . , Am) = - E ^j-T + :,{N - I) J + 25J,. (17) 

i=i ^ + 4 ^ 
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III. GROUND STATE 



Since the parameter Ji depends only on c^, c may take either real or imaginary values. In this section, we discuss 
the ground state properties of different c values. 

(i). c > 5 + 1/2 case. In this case, the coupling between the bulk and the impurity is ferromagnetic. All the A 
modes take real values in the ground state to minimize the energy. Taking the logarithm of (16) we obtain 

r 1 1 

^ = -{^i(A,) + j^m^) - E ^2(A, - AO]}, (18) 



where 0„(A) = 2tan ^(2A/n), 0(A) = 202(s+|c|) (A) — 202(|c|-s)(A) + 02(A) + 0i(A), Ij are some integers. Note above 
we have included the negative modes by putting A^ = — A_j (Aq = 0). Define 

1 1 

Za,(A) = -{^?i(A) + — [0(A)- E ^^2(A-A0]}. (19) 



Then Z^iXj) — Ij/N gives the Bethe ansatz equation (18). For the ground state, Ij must take consecutive numbers 
around zero symmetrically. A density function for the ground state in the thermodynamic limit can be defined as 

PnW - (20) 

with the condition J^^PN^dX — (2M + l)/iV, where A is the cutoff of A modes. As discussed in many earlier 
papersii, the eigenenergy is minimized at A = cx) up to the order 0{N^'^). With this condition, we obtain that 
M = N/2, which gives the self magnetization of the ground state as 

Mg = ^N - M + 2S = 2S. (21) 

Such a result indicates that the impurity moments can not be screened due to the ferromagnetic coupling between 
the bulk and the impurity. 

(ii). 5 < |c| < 5 + 1/2. In this case, > and thus the exchange interaction between the impurity and the bulk 
falls into the antiferromagnetic regime. Two imaginary modes of A at A = i{\c\ — S) can exist in the ground state. 
Note this mode carries energy = — l/[l/4— (|c| — 5)^], which is smaller than those carried by any real modes. The 
real modes thus satisfy the following Bethe ansatz equation 

/' Aj + ^\2N fXj ~ ic + iS ^2 f Xj + ic + iS ^2 fXj — i\c\ + iS + i^2 ^ Xj + i\c\ — iS + is2 



Xj — ^ Xj—ic — iS Xj+ic — iS Xj — i\c\ + iS — i Xj + i\c\ — iS — 



M-2 



Xj — rXi + i 



n n 



(22) 



With the same procedure discussed in case (i), we obtain M — (iV + 2)/2 or the residual magnetization of the ground 
state Mg — 2S □ 1 as expected. This means the impurity moment is partially screened, a similar result to that of the 
Kondo problemD. 

(iii) . < |c| < S. In this case, the system falls also into the regime of antiferromagnetic coupling {Ji > 0). No 
bound state can exist in the ground state. The function Z^^iX) can be defined in a similar way of the case(i) but with 
a different 0(A) = 26*2(5+0) (A) + 26'2(s_c)(A) + 6*2 (A) + 6li(A). With the condition Zn{±oo) = ±{M + l/2)/iV we have 
again M = {N + 2)/2 and the residual magnetization Mg — 2S — 1. 

(iv) . When c takes an imaginary value, Ji is always positive. Suppose c = ib. The Bethe ansatz equation (15) then 
becomes 

+ ^ ■^2N f Xj - b + iS .^2i Xj + b + iS .^2 _ -i-r Xj-rXi+i , 
^X~I> ^X~b^> ^XiTb^' - 4.fjA,-rA,-z- ^''^> 

J 2 r—±l l^j ■' 

For the ground state, all A take real values and their cutoff A still tends to infinity in the thermodynamic limit. As 
discussed in case (i), a similar function Zm{X) can be defined but with a different 0(A) = 26'2s(A — b) + 202s(A + 6) + 
02(A) + 01 (A). Therefore M = [N + 2)/2 for the ground state, which still leaves a residual magnetization of 25* — 1. 

From the above discussion we conclude that the impurities can be screened (partially for S > 1/2) only in the case 
where Ji > 0. The impurity moments can not be screened when Ji < 0, unlike the situation for the Kondo problem 
in a Luttinger liquid predicted by Furusaki and NagaosaS. 
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IV. RESIDUAL ENTROPY AND LOW TEMPERATURE SPECIFIC HEAT 



The-|thermodynamics of the presentjtiodel can be constructed with the standard method proposed by Yapg and 
Yan^y based on the string hypothesisEil. Here we omit the details which can be found in some nice reviewa23. An 
interesting feature of our model is that the residual entropy Sg may take different values depending on the bond 
deformation between the bulk and the impurities: 

i^ii]^, ifici>5+i, 

Sg ^ { ([2| c|]-25)(2-[2|c|] +25y' U "3 + 2 > 1^1 > ' 

21n^4S'2 - [2|c|]2, ifS'>|c|>0 
21n(2S'), if c imaginary, 

where [2|c|] is the maximum integer equal or less than 2\c\. For |c| < 1/2 or imaginary and S = 1/2, Sg = 0, 
which means the impurity spin can be completely screened in the ground state and thus the system flows to a local 
Fermi liquid fixed point at low energy scales. TheJow temperature thermodynamics can be formulated based on the 
local Fermi liquid theory for the Kondo problemEj. Since the spectrum are described by the quantum numbers Ij, 
PNi^j) — 7rZjv(AjXcan be treated as the momenta of the quasi-particles in the Luttinger- Fermi liquid picture for the 
integrable modelaSil. In the thermodynamic limit, the ground state energy can be expressed as 

Eg ^ J eoiX)pNiX)dX + const., (24) 

up to the order of 0{N~'^), where eo(A) = —J/{\^ + 1/4). From the definition of pn{X) in (20), we can rewrite (24) 
as 

^9 = ^ j e{\)p^^\\)d\ + const., (25) 

where e(A) = cos\i(lx) dressed energy function, which can be treated as the energy of the "quasi-particles". Note 

that pn{^) for c imaginary can be solved up to 0{N^^) as 

pjv(A) = po(A) + ^P»(A) + ^P6(A), 

= I XV (26) 
cosh(7rA) 

^"(^^ = 2cosh(vrA) + J 4^cosh(ic.)'^''' 
P. (A) ^ ^ 



cosh(A — h) cosh(A + b) ' 

where po(A), pi{X)/N, and pb{\)/N are the contributions of the bulk, the impurity and the open boundary to the 
density respectively. The density of states at the Fermi surface can be determined in the Fermi liquid picture as 

iSr(n\- ^ dPN^ ^ _ Pjv(A) 

^(°) - 2^^^'^=°° - MA) (^^) 
Therefore, the impurity contribution to the low temperature specific heat reads 

a = J^cosh(7r6)T, (28) 

While for |c| < 1/2, pi{X) can be solved by a similar way as 

^^(^^ " cosh(A - ic) + cosh(A + ic) ■ (^^^ 



Thus the impurity specific heat reads 



= 3^cos(^c). (30) 
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Notice that the Kondo temperature is nothing but the effective Fermi temperature in the local Fermi liquid theory. 
Therefore, we conclude that the Kondo temperature for the cases discussed above is given by 

Tfc = 7rJcos-i(7rc). (31) 

Such a result directly shows the crossover from an exponential law to an algebraic one when c goes from imaginary to 
real as pointed out earlier by Lee and ToneJa. Notice that c = ^Jl — J / Ji and for the Kondo problem in a Hubbard 
chain, the band width of the spinons (proportional to J) is about 4i^/[/. Where 2t is the band width of the fermions, 
U is the onsite Coulomb repulsion. For the strong correlation limit U >> t, J < Ji, the Kondo temperature shows 
a power law of /[UJi). While for the weak correlation limit U « t, J > Ji, the Kondo temperature shows a 
exponential law of /U Ji. 



V. CONCLUDING REMARKS 



In conclusion, we establish an exactly solvable model of Heisenberg chain coupled with two impurity moments. 
This model is relevant to the Kondo problem in a Luttinger liquid. With the algebraic Bethe ansatz method, the 
hamiltonian is exactly diagonalized. It is found that the local moments can be screened (partially for 5 > 1/2) only 
in the case where the coupling between the bulk and the impurities is antiferromagnetic (Ji > 0). Such a result 
strongly contradicts to that of the classical system (Ising model) , where the ground state is a pure Neel state with 
total magnetization of zero (provided N even) for both > and Ji < and arbitrary S. The present result shows 
that the quantum fluctuation plays a crucial rule in the one-dimensional quantum system. The residual entropy of 
the ground state is derived from the thermal Bethe ausatz. It strongly depends on the parameter, c, a similar result 
to that of the Kondo problem in the Kondo pMblemEl. Based on the local Fermi liquid theoryt^ and the Landau- 
Luttinger description for the integrable model^j, we derive the low temperature specific heat for some special cases. 
The Kondo temperatures are exactly carried out for these cases. The exact result does show the crossover of the 
Kondo temperature from an exponential law_to a power law when the parameter c goes from imaginary to real, a 
phenomenon first obtained by Lee and TonerB. 
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